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ABSTRACT 



We use data from the Pan- Andromeda Archaeological Survey (PAndAS) to 
search for evidence of an extended halo component belonging to M33 (the Tri- 
angulum Galaxy). We identify a population of red giant branch (RGB) stars at 
large radii from M33's disk whose connection to the recently discovered extended 
"disk substructure" is ambiguous, and which may represent a "bona-fide" halo 
component. After first correcting for contamination from the Milky Way fore- 
ground population and misidentified background galaxies, we average the radial 
density of RGB candidate stars over circular annuli centered on the galaxy and 
away from the disk substructure. We find evidence of a low-luminosity, centrally 
concentrated component that is everywhere in our data fainter than /ly ~ 33 mag 
arcsec"^. The scale length of this feature is not well constrained by our data, 
but it appears to be of order 

'^exp ^ 20 kpc; there is weak evidence to suggest 
it is not azimuthally symmetric. Inspection of the overall CMD for this region 
that specifically chps out the disk substructure reveals that this residual RGB 
population is consistent with an old population with a photometric metallicity of 
around [Fe/H] ~ -2 dex, but some residual contamination from the disk substruc- 
ture appears to remain. We discuss the likelihood that our findings represent a 
bona-fide halo in M33, rather than extended emission from the disk substructure. 
We interpret our findings in terms of an upper limit to M33's halo that is a few 
percent of its total luminosity, although its actual luminosity is likely much less. 

Subject headings: galaxies: evolution - galaxies: halos - galaxies: individual 
(M33) - galaxies: spiral - Local Group 
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1. Introduction and Background 



A Cold Dark Matter (ACDM) cosmology predicts that larger galaxies are built through 
the hierarchical merging of smaller galaxies. Infalling components are disrupted partially or 
entirely, and part of that material forms the stellar halo of the larger galaxy. Stellar halos 
therefore contain the remnants of past interactions between galaxies, and th eir properties 



can i ndicate the approximate time, size and frequency of past mergers (e.g., 
200^ . 



Pur cell et al. 



We can directly observe only a few relatively nearby halos in any great detail due 
to their faint nature, and only a small number of halos are directly observed through 
their resolved stars outside of the Local Group. Due to their faintness, it is problematic 
to determine whether or not the halos are smooth and/or symmetric. There is likely a 



continuum of sc enarios that we observe 
shells, etc.; e.g., iMartmez-Delgado et al. 



jetwee n newly-accreted objects (creating streams. 



20101 ) and smooth halos, and our interpretation 



will depend on the time since the accretion and the spatial resolution and depth of the 
observations. The l ong dynamical timesc ales for structures outside of the disk implies that 



they are long-lived (IJohnston et al. 



19961) 



Outside of the Local Group, halo detections are extremely challenging as it becomes 
more difficult with increasing distance to distinguish th e halo from other stellar components 



le.g- 



Dalcanton fc Bernstein 



of kinematical data ( [Barker et al 



2002 



de 



one et al. 



2009 



20081 ) and even more so in the absence 



20121). S cattered 



de Jong 



ight and non-stellar pollution 



20081 ). (The surface brightness 



of counts also interfere with halo detections (e.g. 
detection limits generally needed are > 7 magnitudes fainter than the sky where the 
"darkest" skies (20th percentile), at Mauna Kea, are fainter than jjy ^ 21.3.)Q 



^http:/ /www. gemini.edu/sciops/telescopes-and-sites/observing-condition-constraints 
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Searches for halos around distant galaxies began with deep observations around single 
galaxies using surface brightness p hotometry and, as recent studies continue to do, focussed 



on la t e-type edge-on galaxies (e.g. . 



Sackett et al. 



1994 



Shang et al. 



1998 



Zibetti &: Ferguson 



2004 



Tikhonov fc Galazutdinova 2005: Buehler et al. 



20071: Ide Jong et all 120071: ISeth et al 



2007 



Rejkuba et al. 



20091: 



Mouhcine et al. 



2007 



2010 



Radburn-Smith et al 



201ll ). Each 



stellar component is revealed more easily in the cross-section rather than the face-on view. 



for a 


lalo signal ( 


Zibetti et al. 


2004 


2011: 


Zackrisson & Micheva 


2011 





de Jong 



2008 



Bergvall et al. 



together before looking 

0; 



201 



Zackrisson et al. 



20111 ). again highlighting the difficulty of detecting halos 



because of their extreme faintness. 



Halos have also been observed around other types of galaxies, not ius t late-type edge-o n 



galaxies: for example, the Virgo Cluster's central elliptical ga . 



nearby starburst gal axies ( 



elli ptical NGC 3379 

see 



Malin et al. 



Bailin et al. 



Harris et a. 



1983 



Peng et al. 



2002 



2011 



Rys et al 



2011 



Rich et al 



axv. M87 flWeil et al. 



19971 ). 



20121 ). the Leo 



20071). and the giant ellip tical NGC 5128 (Centaurus A; 



Rejkuba et al. 



20111). 



§ [2] provides a literature review of stellar halos in the Local Group. Details of the 
PAndAS observations around M33 are given in § [31 We are ultimately concerned with 
identifying the RGB stars in the M33 halo (if it exists). However, we must exclude the 



re gions associated wit 



m 



1 the extend ed optical substructure surrounding the disk identified 



McConnachie et al 



(12009 . 



20101 ). and we must also correctly account for and subtract 
off the contribution from the foreground Milky Way disk and halo components, and the 
background galaxies misidentified as stars. § H] describes these corrections and exclusions as 
part of the analysis. We discuss our results in § [5l before summarizing in § [61 
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2. Halos of the Local Group Galaxies 



2.1. The Milky Way Galaxy 



The Local Group provides the closest opportunity to study a stellar halo but even the 
Milky Way (MW) is problematic to observe because of the restrictions and biases associated 
with viewing our Galaxy from within - alth ough it has obviously been studied in depth 

and references therein). Current seemingly 



(e.g., see the annual review by 



Helmi 



2008 



contradi ctory evidence means it i s unclear whether th e MW stellar halo is oblate, prolate or 



triaxial f iNewberg fc Yanny 



2006 



Deason et al 



(DM) halo seem to favour triaxiality (iLaw et al. 



20111). although models of the dark matter 



2009 



Law &: Majewski 



2010|). Numerous 



detections of substructure beyond the stellar bulge and disk are anothe r reason that this 



ambiguity remains - substruc ture such as the Sagittarius dwarf galaxy (|lbata et al 



and associated tidal streams fllbata et al. 



ring (llbata et al. 



flMartin et 



Juric et al 



al 



2003 



2008 



Martin et al 



Yanny et al. 



2004ajjb|) and Virgo ( jVivas et al 



2003 



2001bl. 



2002 



Maiewski et al. 



Crane et a. 



2001 



1994h 



20031 ) ■ the Monoceros 



20031). overdensities in Canis Major 



Newberg et 



clouds in the Triangulum- Andror neda region f)Rocha-P i nto et al. 



20071) and t he Hercules-Aqui 



Orph an (iGrillmair 



2006 



Belokurov et al 



a region ( 



2006 



al 



2002 



Xu et al. 



Belokurov et al 



2007bf ) and Cetus Polar (INewberg et al.l 



2006 



2004 



2007ah. and finallv t he 



20091 ) streams. 



There is growing evidence to suggest that the MW halo has a dual halo, with the 



different components a result of their different formation proce s ses (e.g., 



2000 



CaroUo et al. 



2007 



Miceli et al. 



2008 



de Jong et al. 



satellite accretion and in-situ format i on (e.g., 



Zolotov et al 



2009 



OseretaL 



20101 : 



BeU et al. 



McCarthv et al. 



20 



2008 



20121 ). 



Chiba fc Beers 



Beers et al 



20 



Schlaufman et al 



2h, such as 



2009 



2012 



The total (dark plus luminous) mass of the Galaxy within 300 kpc is estimated to be 
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in the range 0.7 < Mmw < 3.4 x 10^^ Mq f lBaiesi Pillastrini 



2009 



Watkins et al 



20101 ) 



The MW stellar halo luminosity, includ i ng all substructure, is est imated to be of order 



^MWMioy 



IO^Lq ( ICarney et al. 



1990 



Bullock &: Johnston 



2005 



an d references t herein) , 



compared to the MW host luminosity, LMw,host,v = S.llgie ^ lO^'^L© ( iSackett 



19971). 



2.2. The Andromeda Galaxy 



Observations of the Andromeda and Triangulum Galaxies (M31 and M33, respectively) 
are free from the problems inherent wit h viewing the MW from within, but are still 



close enough to resolve individual stars ( iMould fc Kristian 



1986 



Crotts 



19861 ) - and 



many groun d-based studies are now also resolving individual stars beyond the Local 



Group (e.g. 



Barker et al. 



2009 



2012 



Bailin et al 



2011 



Tanaka et al 



20111 ). However, 



progress on the study of MSl's full extended stellar halo has only come rel atively recent 
- and the degree to wh i ch it is simi 



2006 



Ibata et al. 



2007 



Koch et al. 



ar to the MW's halo is still uncertain flKalirai et al. 



y 



20081 ). MSl's total (dark plus l uminous) mass is 
estimated to be in the range 9.0 x lO^^Mfr , < M m.^i < 2 x lO^^j^f^ flEvans fc Wilkinson 



2000 



Chapman et al 



2006 



Watkins et al. 



20101 1. The lower limi 



t is determined 



rom 



20061 ). 



the kinematics of RGB stars out to 60 kpc with 99% confidence (j Chapman et al. 
An alternative dynamical mass that uses kinematics of MSl's giant stream rather than 
those of the satellite galaxies, glob ular clusters, plan etary nebulae or RGB stars, gives 



MM3i,R<i25kpc = 7.5l:li X l O^^M^ fllbata et al. 
LM3i,host,v ~ 2.6 X IO^^Lq ( Ivan den Bergh 



20041 ). MSI has a total host luminosity of 



the MW. Many photometri c substructures have been revealed around MSI (llbata et al. 



2001a; 
200% . 



Ferguson et al. 



2002 



Irwin et al. 



2005 



Ibata et al. 



2005 



2007 



MSI ( 


Ibata et al. 


McConnachie et al. 



Irwin et al.l (|2005[ ) fit minor- axis profiles of a de Vaucouleurs law out to a projected 



radius of ~20 kpc, and beyond this a power law (index ~ -2.S) or exponential (scale length 
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14 kpc). 



Ibata et al 



(120071 ) also find a smooth underlying component to which they fit a 



Hernquist profile (scale length ^ 55 kpc), a nd a power law (index ~ -1.91 ± 0.11), similar 



Ibata et al 



estimate that the total luminosi ty of the smooth 



to the MW halo. If symmetric. 

halo is Lmzimio,v ~ lO^i^o, again similar to the MW halo. More recently, 



Courteau et al. 



(I2OIII ) combine ground- and space-based data from several sources and decompose the 
resulting composite luminosity profile to find a halo component described by a power law 
(index ~ -2.5 ± 0.2). M31's halo h as further similarities to that of the MW's in terms of 



metallicity an d velocity dispers ion ( 1 Chapman et al. 



of ~800 stars. 



Chapman et al. 



20061). Using a Keck/DEIMOS sample 



find a non-rotating metal-poor ([Fe/H] ~ -1.4) smooth halo 
between 10 - 70 kpc with no metallicity gradient underlying the metal-rich ([Fe/H] ~ -0.9) 
rotating extended component. Along with a comparable dark matter halo mass, the similar 
metallicities and dispersions are suggestive that the early formation periods of both were 
also similar. 



2.3. The Magellanic Clouds 



The Large Magellanic Cloud (LMC) is the fourth most naassiy e Local Group member . 



It has a total mass of M, 



LMC 



W^Mq fivan der Marel et al. 



2002 



Bekki &: Stanimirovic 



20091). 



20091 ) . The total host luminosity is LiMC,hosty = 3.0 x 10^ Lq (iBekki &: Stanimirovic 
It may also have a stellar halo. A study of RR Lyrae stars, generally associated with old 
and metal-poor stell ar populations, finds that these stars have a large velocity dispersion 



of 53 ± 10 km s"^ f Minniti et al. 



2003 



Alves 



2004j ). A photometric and spectroscopic 



survey has also revealed individual RGB stars enveloping the LMC out to 



and consistent with a de Vaucouleurs profile, suggestive of a classical halo (IMajewski et al. 



arge distances 



20091 ). If the LMC does have a stellar halo it seems somewhat surprising given that another 



tracer of old and metal-poor populations - globular clusters - show no evidence for a halo as 
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they lie within a disk region around the LMC (iFreeman et aLlll983t ISchommer et al.lll992l ) 
This inconsistency could be explained if the GCs are accreted in earlier, mor e gas-rich 



merger events compared to those that populate the halo with individual stars (IBekki 



2003). 



The stellar outskirts of the Small Magellanic Cloud (SMC) have recently been 



Nidever et al. 



20111 ) to reveal a 



studied through the MAgellanic Periphery Survey (MAPS; 
population of nearly azimuthally-symmetric RGB stars out to a radius of ~11 kpc. The 
profile of these stars are well fitted with an exponential profile (scale length ~1 kpc) out 
to ~8 kpc, with a shallower profile beyond (scale length ~7 kpc) - the latter of which the 
authors suggest could be a stellar halo or a population of extratidal stars. 



2.4. The Triangulum Galaxy 



M33, the Triangulum Galaxy, is the thirc 



with a mass close to one tenth that of M31 f lCorbelli fc Salucci 



most massive galaxy in the Local Group, 



20001 measure M33's 

rotation curve out to 16 kpc, and find an implied dark halo mas s of Mm:k?. ^, 5 x 



10 M f7^). M33 has a total host luminosit y of L 



1991) and is more face-on (i = 56° =1= 1°, 



Rubin et al 



1973 



Zaritskv et al 



Athanassoula Sz Beaton 



( Ide Vaucouleurs et al. 



nuclear cluster 



perhaps a bar (jjavadi et al. 



IO^Lq ( Ide Vaucouleurs et al. 



19891 ) than M31 {i = 77°; e.g. 



20061 ). It is classified ClS db SA(s)cd II-III galaxy 
19911). has little or no bulge comp onen t, a UV- and X-ray bright 



Long et al 



981 



Dubus et al 



1999 



and 



Foschini et al 



20041 ). and 



20111 ) , so most of the central light is di s tributed over an 



exponential dis 



McLean &: Liu 



i component (e.g.. ide Vaucouleurs 



1996 



Corbelli et al. 



2008 



1959 



Bothun 



Kormendy et al 



2010 



992 



Minniti et a. 



Javadi et al. 



distance estimates to M33 cover a wide range from 730 ± 168 kpc ( jBrunthaler et al. 



2011 



1993 



I The 



to 964 ± 54 kpc ( iBonanos et al 



20051) 



20061 ). This disagreement appears to arise because of 



the combination of the different techniques used, and also perhaps due to inhomogeneous 



- 11 - 



interstellar extinction in M33. Consistent with iMcConnachie et al.l (120101). we adopt a 



dista nce modulus for M33 of (m - Mq) = 24.54 ± 0.06 (809±24 kpc; 



McConnachie et al. 



20041 ) throughout this paper. This distance is based o n the TRGB method, and is consistent 



Conn et al. 



20121 ). At 809 kpc, 1° corresponds to 14.1 



with our new findings (820t?^ kpc; 
kpc. 

Is there an M33 halo akin to those found in the MW, M31 and possibly the LMC? 
Previous claims of a detection for M33's halo have come from various sources, in studies 
that resolve individual stars (RGB or RR Lyrae stars), or globular clusters, which we briefly 
review here. 



Mould fc Kristian 



( 1l986l ) used the Hale Telescope/PFUEI to observe two fields both 
7 kpc from the centre of each M31 and M33, along their southeast minor axes. By 
comparing the observed giant branches to those for M92 (< [M / H] >= —0.6) and 47 Tuc 



(< [M/H] >= —2.2), they inferred the presence of inner halos. However, 
observed a field with WIYN/S2KB that included the region studied by 



Tiede et al. 



(120041) 

Mould fc KristiaJ 



and found that the peak of the MDF showed a radial varia tion with a gradient c onsistent 



Chandar et al. 



(12002[ ) obtained 



with that of the inner disk region and not of an inner halo. 
WIYN/HYDRA spectra for 107 of M33's star clusters. These clusters, from a sample with 
known integrated HST/WF PC2 colors, were selec t ed to cover the entire age range of M33's 



clusters (6 Myr to >13 Gyr; 



Chandar et al. 



20011 ). 



Chandar et al 



observed a large velocity 



dispersion that, with Monte Carlo simulations, suggested that old (> 1 Gyr) clusters could 



be split into two c omponents which 



a halo. Similarly, 



Sarajedini et al. 



they associate with a disk population and the other with 



(120061 ) observed 64 RR Lyrae variable stars (type RRab) 
using HST/ACS to have a double peak in periods, again suggesting two subpopulations: a 
disk and a halo component. RR Lyrae stars are only observed in populations older than 10 
Gyr, and therefore stars in these populations should be at least as old as the RR Lyraes. 
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Fu rther hints of M33's halo also come from the following sources. [Barker et al 
(12007a 3) mspected three H ST/ACS southeastern field s ~20'-30' (~4. 7-7.1 kpc, assuming a 



distance to M33 of 809 kpc; 



McConnachie et al 



20041 ) from M33's nucleus. Mixed stellar 



populations were revealed in the CMDs, with an age range from < 100 Myr to a few 
Gyr. The authors compared synthetic populations with the observed CMDs, and found 
that the mean age increased with radius from ~ 6 to 8 Gyr, and the mean metallicity 
decreased from ~-0.7 to -0.9 dex. They concluded that while the fields ar e dom inated 



Cioni et al. 



(120081) used 



by a disk population, a halo component may also be present. 
UKIRT/WFCAM near-infrared observations of a 1.8 degree^ region centred on M33 to look 
at the ratio of C- to M-type AGB stars. They also found metallicity and age gradients 
such that the outer regions were more metal po or a nd a few Gyr ol der than the central 



regions, in agreement with 



Barker et al. 



fl2007a 



b|). iFergusonI (120071 ) reviews results for 



M31 and M33, and notes that while a halo-like component with a power-law structu re wa s 



proving elusive, the RGB narrows and becomes more metal poor beyond ~10 kpc. 



Teigl 



( 120081 ) used RGB and AGB star counts along M33's minor axis, and also observed a break 
in the surface brightness profile at 11 kpc. The pr o file at this region appeared to change 



from an exponential to a power-law. 



Barker et al. 



(1201 Ih use HST/ACS to observe two 



fields 9.1 kpc and 11.6 kpc along M33's north major axis. They find that the outer field is 
old (7 ± 2 Gyr), moderately metal poor (mean [M/H] ~ -0.8 ± 0.3), and c ontained ^30 



times less stellar mass than the inner field. One of the interpretations that iBarker et al 



discuss is that the outer field is a transition zone from the outer disk to another structural 
component. 



Grossi et al. 



( 1201 ll ) use Subaru/Suprime-Cam data with seven fields 10 < r < 
30 kpc from the centre of M33 in the NW and SE. An exponential scale length of ~7 kpc is 
found for both regions, and these authors favour that this component is an extended disk 
rather than a halo. 

The previous section highlights that unambiguous detections of the various galactic 
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components - even for one of our closest neighbouring galaxies - are still extremely difficult. 
While studies of star clusters and RR Lyrae stars show evidence for an old halo, studies of 
individual stars in other subpopulations have only been marginally conclusive. Part of the 
problem was the limited coverage and/or depth. 



McConnachie et al. 



( 120061 ) undertook a spectroscopic survey of RGB stars using 



Keck/DEIMOS. Ra dial velocity distribu tions of these stars were best-fit by three Gaussian 



components, which 



McConnachie et al. 



interpreted as contributions from a halo, a disk. 



and a component offset from the disk (whi ch they suggested could have been a stellar 



st ream or another stella r halo component) 



Ibata et al 



of 



Ferguson et al. 



( 120071 ) extended the observations 



(120071 ) with CFHT /MegaCam along the southeastern corner of M31's 



halo out to M33's centre. 



Ibata et al. 



clearly saw the classical disk of M33, and in addition 
revealed an extended component. They fit a profile to the data between 1 and 4 degrees (the 
edge of the disk, and the point just before where the profile starts rising again, respectively) 
and found the exponential scale length to be 18 ± 1 kpc, or 55 ± 2 kpc using a projected 



Hernquis t model. T 



although 



Ibata et al. 



lese scale lengths were surprisingly as big as they found for M31, 
cautioned that without a full panoramic view it was not possible to 



determine whether or not this feature was a "bona fide" halo. 

Direct and unambiguous evidence for a stellar halo around M33 remains elusive but is 
th e aim of this paper. Any such com ponent must be quite faint. We extend the work begun 



by 



McConnachie et al 



( 12009 



20101 ) as part of the Pan- Andromeda A rchaeological Survey 



(PAndAS), which itself 



juilt on previous 



and CFHT/MegaCam (jibata et al. 



surveys with the INT/WFC ([Ferguson et al 



20071 ) 



20071 ). Optical observations prior to PAndAS suggested 



M33' s disk had an undisturbed appearance - a view that persisted until relatively recently 



(e.g.. 



Sarajedini 



2007 



Ferguson et al 



20071 1. This implied that the disk had not been 



tidally disrupted by either the MW or M31, and was seemingly discrepant when compared 
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to the radio detecti on of a warped gaseous disk component (e.g., iRogstad et al.lll976 : 



Putman et al. 



Roestad et al. 


1976 




with the 


Ferguson et aL 



data set was depth and that the observations could only rule out the presence of substructure 
or a halo with surface brightnesses brighter than fiy ~31 mag arcsec"^. 

PAndAS observations covering the area around M33 with unprecedented combination 
of depth and coverage have revealed a vast low surface brightness stellar substructure. 
The S-shaped optical warp of this substructure is generally aligned with the HI warp, and 
therefore resolves the previous discrepancy. It now s eems most likely tha t this warp was 
the feature that was previously partially detected by 



McConnachie et al. 



castin g doubt on the previous interpretations from this kinematic study) and 



fl2006r) ftherebv 



Ibata et al. 



(120071 ). Although the nature of the substructure is still being investigated, the favoured 



interpretation for its origin is that it is a disruption of the dis 



interaction with M31 as M33 orbits M31 (iMcConnachie et al 



m prep 



i that was caused by a tida l 



2009 



2010 



Dubinski et al 



). Preliminary models can reproduce the s hape of the extended d isk substructure, 
and also satisfy M33's proper motion constraints fiBrunthaler et al.li2005i'). Spectrosco pic 



observations may provide further clues (see upcoming paper by 



Trethewey et al. 



20121). 



The depth of the PAndAS data allows us to test whether or not a genuine halo 
component is observable in addition to the disk-like warp. Given that the warped extended 
disk substructure is extremely faint and had eluded detection for so long, it would be 
reasonable to expect any underlying component of a halo to also be extremely faint. Indeed, 
a faint extended stellar component is hinted at beyond the e xtended disk substructure , and 



one possibility put forward is that this is a halo component (IMcConnachie et al. 



20101 ). 



Since M33 is r elatively low mass, would we expect it to have a detectable halo? 



Purcell et al. 



( 120071 ) predict that a galaxy with total mass M ~ 10 will, on average 



have a halo that contributes <1% of the total luminosity from the galaxy. Therefore, the 
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expected total halo luminosity of M33 could be as low as Lmssm^o 

< WLq. Their halo 

estimates make no distinction between the smooth component and the substructure. Here, 
we define the "halo" as the component or stellar population in the outer regions of M33 
th at is not clearly associated with the disk or the extended disk substructure identified 



by 



McConnachie et al. 
^ 



torn. N 



halos, similarly to iPurcell et al. 



ote t hat we do not distinguish between smooth or lumpy 



(120071 ). 



3. Observations, Data Reduction and Calibration 



We use data 



McConnachie et al. 



Tom t he Pan- Andromeda Archaeological Survey (PAndAS; 
2OO9I ) to observe 48 degree^ around M3 3 with CFHT/Meg aCam 



Martin et al. 



20091 ). The data 



out to a projected radius of 50 kpc (cf. rM33,viriai = 152 kpc, 
has limiting magnitudes for point-source detections of g' ~ 25.5, i' ~ 24.5 (AB magnitudes 
on the SDSS scale) at S/N = 10 in subarcsecond seeing. MegaCam is composed of 36 
individual CCDs, has a 0.96 x 0.94 degree field of view and a resolution of 0.187" pixel"^. 

Figure [1] shows each location of the MegaCam ~square- degree images around M33 in 
a tangent-plane projection. To be explicit, for our analysis we only consider the data for 
MegaCam images within the annulus with radius 3.75 degrees centred on M33 as shown 
in Figure [1] The MegaCam subexposures are dithered in order to cover the small gaps, 
but not the large gapj^. In the large gaps there may be fewer detections due to shallower 
depths. However, the area of the chip gaps is very small compared to the overall MegaCam 
field and can be neglected for the purposes of this analysis. 

The prefixes of the image labels in Figure [1] represent the timeline of the observations: 
The central field, m33c, was observed primarily in the observing semester 2004B and 



See http : / / www. cfht . hawaii . edu / Instruments / Imaging / Megacam / specsinf ormation . html 
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retrieved from the CFHT archive, with some data from 2003B. All other fields with prefix 
m were observed in 2008B. Fields with prefix nb were o bserved in 2009B. D ue to a failure of 
CCD 4 in the 2003B observing semester, the data from 



Ibata et al 



(120071 ) which extended 



the southeastern section of M31's halo in a line to the centre of M33 was replaced with 
data from 20 10B (prefix tb). The ellipse in Figure [T] marks the fiB ~ 25 mag arcsec"^ 



( iNilson 



19731 ) contour of M33's disk, and the solid-line cross represents themaj or and minor 



axes of M33, with the major axis inclined 23 degrees to North ( INilson 



19731 ) ■ The dashed 



concentric circles represent radii at r = 1, 2, 3 and 3.75 degrees (14.1, 28.2, 42.4 and 53.0 
kpc, respectively) from the centre of M33. The data within the annuli that they delineate 
will be used in the analysis that follows. M33 is approximately 31 degrees below the central 
axis of the Milky Way disk, M33(i^f,) = (133.61, -31.33) degrees, compared with M31 which 
is about 21 degrees below, M31(/,fe) = (121.18, -21.57) degrees. The three dashed lines in 
Figure [Hare lines of equal Galactic latitude (6 = -35.3, -31.3 and -27.3 degrees). 

Pre-processing and reduction were undertaken with Elixijf] by the CFHT team, and by 



the Ca mbridge Astronomica. 



i mage s 



(12009 



Irwin fc 



20101 1 and 



2001 



Survey Unit (CASU) through a pipeline adapted for MegaCa m 



), respectiv ely. The reader is referred to 



McConnachie et al. 



Cockcroft et al. 



(120111 ) for more details. 



4. Analysis 

Taking advantage of the wide coverage of the PAndAS data, we deliberately seek 
direct evidence for M33's stellar halo in this data, and expect it to be extremely faint, 
centrally-concentrated, and detectable via RGB stars. However, this low-luminosity 
component will be mixed with stars from the M33 disk, M33 extended disk substructure 



^http:/ /www. cfht.hawaii.edu/Instruments/Imaging/MegaPrime/dataprocessing. html 
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surrounding the disk, and the MW foreground (both its thick disk and halo), in addition 
to background galaxies misidentified as stars. Our technique involves statistically removing 
the MW foreground stars and background galaxies, excluding the regions identified as 
belonging to M33's extended disk substructure, and seeing what signal remains. 

Figure [2] shows the colour-magnitude (Hess) diagrams for the data in the annuli 
in Figured] We note that Figure |2] contains more than 1.4 million objects that were 
identified as robust stellar candidates in both go and io through the CASU pipeline's object 
morphological classification. Magnitudes are de-reddened source by source using values of 



E{B - V) in the range 0. 034 < E(B - 
io = i- 2.086 E{B - V) Jschlegel et al. 



V) < 0.130, with go = g - 3.793 E{B - V) and 



19981 ). The data is binned in 0.025 x 0.025 mag 



bins and is shown with a logarithmic scale for the number counts of stars. As mentioned 
previously, we want to identify M33 RGB stars but first we need to estimate the level of 
contamination. To examine the MW foreground contamination, we look at the two sources 
of contribution from MW stars easily identifiable in the CMDs. The MW halo turn-off 
stars are seen as a thin band on the left of the CMDs, and we use a region defined as 0.1 
< {9 ^ < 0.6, 19 < io < 22 to measure their relative numbers in each zone. The red MW 
disk dwarfs are seen as a broader band on the right, and we identify them in the region 
1.5 < {g — i)o < 3, 17 < in < 20. Both of the regions for the MW disk and halo stars are 



consistent with 



McConnachie et al. 



020 1 01) • Finally, M33 RGB stars are selected by the 
colour- magnitude locus where we would expect to find RGB stars . This locus is def i ned 
using isochrones from the Dartmouth Stellar Evolution D atabase ( 



which are transformed to the CFHT photometric system (IMcConnachie et al. 



Dotter et a 



2007 



20081) 



20101). These 



isochrones are between the 12 Gyr [a/Fe] = 0.0 isochrones, shifted to the M33 distance 
modulus, with metallicities of -2.5 dex < [Fe/H] < -1 dex. This is a necessarily broad cut 
to allow for the possible range of metallicities that may be present in M33's halo, which we 
expect to be predominantly metal poor. Note that metal-rich stars may also be present. 
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but will likely contribute a small amount to the overall halo component, while increasing 
dramatically the contamination from foreground stars that occupy a similar locus in the 
CMP . We could expec t some a-enhancement in the M33 halo, as we see in the MW halo 



(e.g. 



Venn et al. 



2004 ). but since there is no evidence to suggest this we adopt [a/Fe] = 
0.0 for simplicity. We also note that the isochrones are being used to help define a locus 
in the CMD, and an absolute interpretation of the implied metallicities is not intended 
(for example, there will also be age degeneracies). A magnitude limit of 21.0 < iq < 24.0 
is also imposed on the RGB candidate stars, with the lower limit ensuring a high level of 
completeness while excluding the majority of bright background galaxies mis-identified 
as stars (which becomes a major source of contamination at faint magnitudes; io ~ 25, 
^ "^0 ^ 1). We test the effect of raising the faint limit to io < 23.5 in Section 14^31 

The four panels in Figure |2] correspond to annuli with the radii between r = 0-1, 1-2, 
2-3 and 3-3.75 degrees. We use the latter annulus to estimate the spatial variation in the 
MW foreground since any M33 halo component, if present, is likely to be very weak. The 
number of stars in each annulus, and the number of stars within each of the three selection 
regions, are shown in Table [H 



4.1. Extended Disk Substructure 



Figure E] is a revised version of Figure 13 in 



McConnachie et al 



f l2010h . using the 



ne w data in images tb 62-tb66 (see Figure [T]) . The map was created in an identical way 



to 



McConnachie et al. 



see their Section 3.2.2 for details). Figure |3] shows the density 
contours of candidate RGB stars, and uses a slightly narrower metallicity cut of -2.0 dex 
< [Fe/H] < -1.0 dex than the cut we impose on the candidate RGB stars in the CMDs. 
This narrower cut is used simply because this is the metallicity range in which the extended 
disk substructure component is strongest. There is hardly any contribution to the extended 
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disk substructure from stars with metallicity between -2.5 dex < [Fe/H\ < -2.0 dex; 
however, we would not necessarily expect this to be true of M33's halo RGB stars. The 
single grey contour represents la above background, or an estimated surface brightness 
limit of iiv = 33.0 mag arcsec"^. The other (black) contours are 2, 5, 8 and 12cr above the 
background (/xy = 32.5, 31.7, 31.2, and 30.6 mag arcsec"^, respectively). 

Figure H] shows the contributions to the total radial profile from the regions defined 
both within and excluding the la contour shown in Figure [3l The profiles for the extended 
disk substructure and non-substructure regions are normalized using the total annulus area. 
The non-substructure regions are seen to start dominating the profile for r > 2 degrees. 
We exclude data within the la contour when probing for the stellar halo. When we excise 
the extended disk substructure area denoted by the la contour, note that we cannot probe 
radii smaller than r < 1 degree. 



4.2. Foreground and Background Contamination 

We have identified candidate stars for the M33 RGB, and MW disk and halo 
populations, and we have identified the regions associated with the extended disk 
substructure surrounding the disk. We now test the populations for variations in the spatial 
distributions. Figure [5] shows smoothed non-excised maps of the spatial distribution for 
the background galaxies and each of the three populations identified in Figure [2] (i.e., the 
MW disk, MW halo and the M33 RGB candidate stars). We identify background galaxies 
morphologically using the CASU pipeline, and those shown in Figure |5] have had broad 
colour and magnitude cuts applied (17 < < 23.5, 17 < < 23.5, and -2 < {g — i)o < 4). 

The data is binned into 18 x 18 arcsecond cells. The galaxy, disk and halo maps 
are smoothed once, and the RGB map is smoothed three times, all with a boxcar size 
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of 13 X 13 cells (or 



McConnachie et al. 



equiv alently 3.9 x 3.9 arcminutes; exactly four times smaller than in 



2010l ). The RGB map is smoothed three times to better highlight the 



faint extended disk substructure surrounding the disk. 

The galaxy, MW disk and MW halo maps clearly show no significant global features, 
although the centre of M33 is apparent due to the crowded nature in this region where 
the automated object morphological classification is less successful. Apparent holes in the 
data are caused by bright foreground stars preventing detection of faint objects in their 
surroundings. The galaxies misidentified as stars in our sample are expected to have a 
similar distribution to the galaxies shown in the galaxy map. The RGB map shows the 
extended substructure surrounding M33's disk, and Andromeda II to the north-west. 

Now we excise the regions associated with the extended disk substructure surrounding 
the disk and investigate the variations of the MW disk and MW halo populations in 
different regions on the CMD within the 3-3.75 degree annulus (in which we expect little 
contribution from bona-fide M33 stars). 

Figure [6] shows the variation of these three populations with respect to the azimuthal 
(left-hand column) and Galactic latitudinal (right-hand column) distributions. All panels 
show the variation between 3 < r < 3.75 degrees, with extended disk substructure regions 
excised. Each of the three rows shows the density variation of M33 RGB, MW disk and 
MW halo candidate stars. The density of disk stars increases towards the disk, as does 
the density of the stars in the MW-halo selection region but with a smaller amplitude. 
As we do not expect the halo stars to vary in latitude in this manner, this suggests some 
cross-contamination with the thick disk stars. 

Within the RGB selection shown in Figure IHl there is little variation in the annulus 
at large radii. Indeed, the best fit weighted least-squares fit in both RGB panels is 
consistent with a slope of zero. As such, we conclude that there is no reason to adopt a 
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spatially- varying foreground for our analysis, and instead use a constant background, Tjbg. 



Having determined the extended disk substructure area to avoid, we produce 
substructure-excised radial profiles. As previously stated, we expect M33's halo to be 
extremely faint and centrally-concentrated so we bin the data in annuli centred on M33, 
where we require a certain signal-to-noise ratio for the bins in each profile. 

Figure [7] shows radial density profiles of the RGB stars after excising extended disk 
substructure regions. The small vertical radial density error bars are calculated using 
-y/n/area as the error on the mean of the star counts in each stellar population. The 
horizontal error bars indicate the width of the annulus. The size of the annulus was allowed 
to vary until the signal-to-noise reached the required value (where the "noise" is the radial 
density uncertainty). Each bin in Figure [7] has a signal-to-noise (S/N) cut of 25. We also 
use different S/N cuts, but later show that the results are statistically the same. The larger 
error bars shown in Figure [7] show the variation due to residual substructure. These latter 
errors were measured as the standard deviation of number counts between azimuthal bins 
(36 degrees in width) around a given radial annulus. 

In all the radial profiles, we see evidence for a low-luminosity and centrally concentrated 
profile in M33's RGB stars, which is beyond the extended disk substructure surrounding 
the disk, and has not previously been seen. For illustrative purpose only, as this component 
is so faint and the error bars are large, we use a Levenberg-Marquardt least squares method 
to fit the following exponential model, as shown by the curved lines in Figure [3 



4.3. 



Radial Profile 




(1) 
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The data points are overlaid with the best fit, as shown by the curved dashed hne. The 
horizontal dashed lines show the background level estimated by the fit. We also show the 
background-subtracted fit with the solid curved lines at the bottom of each panel, where 
the use of a constant background is justified in the previous section. Table [2] shows the 
parameters associated with each of the fits at different S/N cuts, including the S/N = 25 
cut shown in Figure [71 As previously mentioned, although the parameters vary slightly for 
each different S/N cut used, they are statistically the same. 

We also test the effect of raising the faint limit of the RGB selection criteria to io = 
23.5 from = 24.0 magnitudes. When we make the brighter magnitude cut we exclude 
more contamination - as we would expect - but we include proportionally less signal. 
The form of the radial profile is essentially the same, although less defined. We therefore 
continue the analysis using the io = 24.0 cut. 



We approximate an equivalent surface brigh tness scale by using 



';he conversion between 



McConnachie et al. 



torn (specifically. 



stars counts and surface brightness described in 
for riRnP! < 350 sta r s degr ee"^ from their Figure 15). For details of this conversion, see 



McConnachie et al. 



(120101). but n ote that the conversion is only an approximation. The 



RGB stars in iMcConnachie et al. 



are selected using -2 < [Fe/H] < -1 dex, and io < 23.5 
magnitudes, whereas here we use -2.5 < [Fe/H] < -1 dex, and io < 24.0 magnitudes. There 
are also large systematic uncertainties inherent in the technique. 

We estimate the luminosity of this component by first simply summing the total 
number of stars contributing to the profile in the radial range for which we have data (i.e., 
0.88 < r < 3.75 degrees), and using the conversion as above. Assuming Poisson statistics, 
we obtain 765 ± 95 stars (assuming a background of 355 stars degree"^, and without 
propagating the uncertainty in the background), corresponding to a luminosity of L = 2.4 
± 0.4 X IO^Lq. Note that this initial estimate is independent of any assumptions we could 
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make about the profile of tlie component. 

To calculate the luminosity extrapolated to the center, however, we assume a spherically 
symmetric smooth profile that is described by the exponential fit with a scale length of 
1.5 degrees (or 21.1 kpc). We calculate the fraction of the integral of the exponential 
fit between 0.88 < r < 3.75 degrees compared to < r < 0.88 degrees, allowing us to 
calculate the luminosity between < r < 3.75 degrees. This simple extrapolation yields 
an estimate of L = 3.8 ± 0.5 x 1O^L0. (If we use a similar technique to extrapolate under 
the whole exponential curve, i.e., out to infintiy, we obtain L = 4.1 ± 0.5 x 1O^L0.) We 
note that the large uncertainties on the exponential profile fit to our data and the unknown 
intrinsic profile of this component, make these extrapolated estimates highly uncertain. As 
noted above, we also do not include the error on the background. The effect of including 
this is seen in Figure [HI as we integrate out to larger radii, the relative luminosity error 
estimates increase. At larger radii, there are fewer candidate RGB stars but a relatively 
larger contribution from the background. 



Discussion and Conclusions 



It was expected that any stellar halo signal around M33 would be at least as 
faint as the recently discovered extended optical disk substructure surrounding the disk 



(iMcConnachie et al 



2009 



20101 ). Hints of a radial fall off beyon d the extent of the 



exten ded disk substructure suggested a tentative halo detection (IMcConnachie et al. 



20101 ). We followed up this possibility in our present study, using higher spatial resolution 
maps, excising any contribution from the extended disk substructure, and subtracting off 
contamination from foreground and background sources, so that we are more able to cleanly 



resolve and identify any remaining signal. 
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We detect a radial density drop off that we interpret as an upper limit of the 
M33 candidate stellar halo. The signal is extremely faint, but seems robust to various 
signal-to-noise cuts; as previously noted, we observe only 765 ± 95 excess stars between 
0.88 < r < 3.75 degrees. 

Are we justified in claiming this extra component is a halo? We have azimuthally 
averaged annuli centred on M33 to find a low-luminosity and centrally concentrated profile. 
The top panel of Figure [9] shows the azimuthal distribution of the RGB candidate star 
density within 3 degrees having excised the extended disk substructure regions. We see 
contamination of the MW foreground stars does not appear to affect the density variation 
of these RGB candidate stars (i.e., we do not see a reflection of the density profiles shown 
in Figure |6] for the MW disk candidate stars). If our extra component was actually residual 
low- level emission from the already known extended disk substructure we would expect to 
see this reflected in this plot, with overdensities around the regions associated with the 
tips of the S-shaped warp, indicated by the two arrows in the top panel. The azimuthal 
distribution is fairly fiat, but we note that overdensities are apparent near the warp's tips 
suggesting some contamination from the extended disk substructure. We further split 
the data into two annuli, 1-2 and 2-3 degrees, but we do not see evidence that the RGB 
candidates show any major differences in their azimuthal distribution from one another in 
either annulus. 

Further constraining this newly discovered component, we show in Figure [10] the CMD 
for all objects with r < 3 degrees, except for those within the extended disk substructure's 
Icr contours in Figure |3] (again, this imposes a minimum radius of ~1 degree). We note 
that the RGB stars that we aim to detect are just visible to the eye on the left-hand plot. 
Again, we see the extreme relative faintness of this component. On the middle panel we 
overlay an [Fe/H] = -2 dex isochrone to this feature. As expected if this component is a 
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halo, this crude measurement indicates that it is relatively n aetal-poor. The extendec 



disk 



2010h . We 



substructure metallicity for comparison is [Fe/H] = -1.6 dex (iMcConnachie et al.l 
show the CMD for the extended disk substructure in the right-hand panel, again overlaying 
an [Fe/H] = -2 dex isochrone for comparison. We can see that the candidate halo lies on 
the metal-poor side of the extended disk substructure RGB. 

With this component that we identify as M33's candidate halo, it is appropriate to ask 
- even with such poor signal-to-noise - if we see any azimuthal asymmetry. To test for this, 
we split the data into the four quadrants split by the major and minor axes, e.g., as shown 
in Figure [31 The resulting radial profiles for each quadrant are shown in Figure [TT], and the 
associated CMDs are shown in Figure [121 It appears as if the east and south quadrants 
have the steepest radial declines, whereas the north and west quadrants are flatter. In other 
words, there may be a variation in the radial profile of the candidate halo when split along 
the major axis. Further interpretation of this possible asymmetry must wait for higher 
quality and deeper data. 

In summary, we have a weak detection that is not clearly indistinguishable in either 
azimuthal distribution or metallicity from the extended disk substructure component. If 
this is a halo, we use the estimates so far obtained to place upper limits on the luminosity. 
We note that with the data at hand we must leave open that it could be another component, 
such as a very extended thick disk. 

The location of the extended disk substructure was the most important knowledge prior 
to beginning this study, in a similar way that the spectroscopic k nowledge of the metal -rich 



component led to the discovery of the metal-poor halo in M31 (I Chapman et al. 



Kalirai et al. 



20061 : 



20061 ). If we directly compare MSl's halo with M33's candi date halo we find 



that apart from the obvious difference in luminosity {LM3i,haio,v ~ 10 i^e 



05 



Ibata et al 



2007 



LM33,haio,v = 4.1 ± 0.5 X IO^Lq.), expectcd because of the mass difference between the two 
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galaxies, it is unclear if the exponential scale lengths are significantly different; we estimate 
M33's scale length to be 21 ± 18 kpc, similar to that found for M31 



20051 ). In light of the discovery here, the spectroscopic work by 



14 kpc; 



Irwin et al. 



McConnachie et al 



mom 



nee ds to be reyisited so that a comparison of M31 and M33's halo 



see 



Trethewey 



2011 



, and a forthcoming paper by 



Trethewey et al. 



metallicity can be made 



20121). 



As mentioned in Section 12. 4[ the favoured interpretation to explai n M33's extended 



disk s ubstr ucture surrounding the disk is a tidal interaction with M31 (iMcConnachie et al 



2009 



20101 ). It is extremely likely that this interaction also affected M33's halo: at least 
altering if not stripping it, with some of M33's halo then being accreted onto M31. The 
halo could also extend beyond the point to which the PAndAS data set is able to measure 



it 
of 



3.75 degrees, or ^ 5 de grees to the north-east.) In the supplementary material movie 



McConnachie et al 



(120091), the end of the modelled interaction also includes the stellar 
halo. Though a significant amount of halo material is stripped from M33, appearing to 
extend beyond the virial radius to form a low-luminosity bridge between M33 and M31, 
most of the halo appears to remain bound to M33. The material stripped from M33's 
halo also extends well beyond the area observed in PAndAS. More of the remaining bound 
material appears on the south-east side (away from M31) than in the north-west (closest 
to M31). Our observations appear to broadly agree with this model, as we see more of a 
gradient in the radial profile in the south quadrant. 

The most probable scenario(s) for how M33's candidate halo was built could be quite 
different from Milky Way and M31 because it is approximately ten times less massive than 
either. Unlike M31, M33 has no bulge, a warped extended disk substructure, and is likely 
interacting with a much more massive neighbour. Studies of the M33 outer halo clusters 
also suggest that with the low GC surface density (Sgc,M33 ~ 0.14 deg^^) compared to 
M31 (SccA/si ~ 0.8 deg~^), M33 either had a much calmer accretion history than M31 or 
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t hat some of the outer halo clusters could have been tidally stripped by M31 ( iHuxor et al. 



2009 



Cockcroft et al.ll201l[ ). The latter idea obviously supports the favoured interpretation 



that could explain the warped extended disk substructure component. 



We now compare our results with a model that predicts the size of M33's stellar 



halo. 



Purcell et al. 



( 120071 ) use an analytic model with empirical constraints from 2; ~ 
observations to predict the fraction of stellar halo mass compared to the total luminous 
mass. They define the diffuse stellar mass fraction as fmi = Mf -^^/Mf*"'; note that 
they use mass rather than luminosities, to avoid uncertainties involved with luminosity 
evolution. There is no distinction made between the substructure in the halo, or the smooth 
diffuse halo that might underlie the substructure. Their predictions cover a range of host 
galaxy's dark matter halo masses, from small late-type galaxies to large galaxy clusters 
(~ 10^^ to ~ 10^^ ^0)- The stellar material is assumed to be able to become part of the 
diffuse stellar halo when its dark matter subhalo has become significantly stripped. The 
dark matter subhalo is considered disrupted when its maximum circular velocity falls below 
a critical value - which is set by considering the em pirical constraints. For DM halos of 



mass ~ 10"Mp 







M; 



M33) 



Corbelli fc Salucci 



20001 ). the stellar halo luminosity fraction is 



expected to be < 1.0% (thus, Lgteiiarhaio < IO^Lq). A galaxy's mass-to-light ratio (over the 
entire halo out to the virial radius) varies as a function of DM halo mass, and this drives 
the fraction of stellar halo material; small galaxies are expected to accrete material from 
dwarf galaxies, which have high mass-to-light ratios and therefo re share little luminous 



material with their host galaxy's stellar halo (jPurcell et al. 



20071 ). Even if they share all of 



their material, the contribution is not large. This picture seems to broadly agree with our 
upper limit estimate of MSB's extremely faint candidate stellar halo, LMssAaioy = 4.1 ± 0.5 
X lO^Le (0.4% < Lm33 ,host,V 0.5%). 



Figure [13] plots the fraction of halo luminosity compared to the host galaxy luminosity. 
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against the host galaxy mass for sev eral galaxies including M33. The lines represent 



equation 7 from 



Purcell et al. 



( 120071 ) for the model of the intrahalo light, with different 
values for the parameters rie// and fa- Here, ^g//, expected to be of order unity, represents 
the effective number of satellites with mass Msat = Mhost/'^'O; fd represents the total stellar 
mass fraction a satellite contributes to its host galaxy halo. We note that if we swap the 
values for fa and Uejf, the lines would vary in the same way. 

The two estimates for M33's candidate halo luminosity fraction are for the directly 
observed estimate (0.88 < r < 3.75 degrees), and the implied, extrapolated estimate 
(r < 10.64 degrees). The values for the MW and M31 are taken from Sections 12.11 and 



(Barker et al. 


2012 




Fraternali et al. 


2002) 



20121). as it has a similar total stellar ma ss to M33 (9.4 ± 0.7 x IO^^Mq; 



20021 ). It is at a distance of 3 .1 Mpc ( iFreedman et al. 



inclination of 63 degrees (IFraternali et al. 



20011), has an 



2OO2I ). and is the brightest member of a loose 



galaxy group and is therefore considered much more isolated than M33 (the closest large 



galaxy is M81, w^ 



Barker et al. 



20121). 



lich is four times further from NGC 2403 than M33 is from M31; 

use Subaru/Suprime-Cam to obtain images 39 x 48 kpc 



Barker et al. 



around the centre of NGC 2403, and see an extended component which could be disk 
structure or a halo. Extrapolating out to 50 kpc they find that the haloes contain ~ 1-7 % 
of the total V-band luminosity, or L2403,haio,v ~ 1-7 x IO'^Lq, depending on whether or not 
an exponential or Hernquist profile is used (if they extrapolate out to 100 kpc the estimate 
does not significantly change). 

The valu e s of L hninl Lhn.'ifnninT.y for the MW and M31 are close to those of the models by 



Purcell et al. 



( I2OO7I ). However, for the less massive galaxy NGC 2403 we find the models 
seem to underestimate the contribution of halo light from these smaller galaxies - and may 
also do the same for M33, although with such a weak signal as we detect here there are 
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large uncertainties. If we further include M33's extended disk substructure {L substructure ~ 
IO^Lq) in the "halo" term, we see that M33 lies even further from the model lines. How do 
we interpret this information? A value of fd = 1 implies that any satellite galaxy has been 
completely destroyed and contributed all of its material to the halo. The MW and M31 
data seem to favour a value of fd = 1, which is of course inconsistent with observations 



'e.g., M31's latest tally is up to 29; iRichardson et al. 



2011 



Slater et al. 



2011 



Bell et al. 



20111 ). If the M33 candidate halo fraction is closer to the upper bounds, then it also appears 



that the models underestimate the halo fraction for lower-luminosity galaxies. 



6. Summary 

We use Pan-Andromeda Archaeological Survey (PAndAS) data to identify RGB 
candidate stars in the regions unrelated to the disk and extended disk substructure 
surrounding the disk. Contamination from both Milky Way foreground stars and 
misidentified background galaxies is subtracted. We reveal a new component centred on 
M33 that has a low luminosity. With such a weak signal, measurements are not well 
constrained by our data. However, it appears that this component has an exponential 
scale length is of order 

"^exp ~ 20 kpc, a photometric metallicity of around [Fe/H] ~ -2 
dex, a luminosity range of a few percent of M33's total host luminosity, and is azimuthally 
asymmetric. More observations and deeper photometry are required to better determine 
the detailed structure of the stellar populations. 

If this feature is truly a halo, it provides support that stellar halos are a ubiquitous 
component of all galaxies, built through the hierarchical merging predicted in ACDM 
cosmology. 
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Table 1. The number of stellar objects located in each annulus shown in Figured! and in 
each region shown in Figure O The sixth column includes all points out to a radius of r < 

3.75 degrees. 







Annuli (degrees) 




Region 


0-1 


1-2 


2-3 


3-3.75 


All 


M33 RGB 


83394 


7597 


7235 


3472 


101698 


MW Disk 


9064 


15431 


26326 


13716 


64537 


MW Halo 


9894 


4309 


6784 


3479 


24466 


Total within annulus 


595878 


163178 


243949 


148885 


1151890 



Table 2. The fit parameters for the exponential model fits in Section H73] shown in Figure 
[71 Each column shows the quantity first in measured units, then in physical units in 
parentheses. Units for each column are shown as footnotes. The uncertainties on each 
parameter are shown under the columns labelled with A. 



S/N So " ASo ro Aro T^tg AJ^tg " 



15 233 (3.7) 189 (3.0) 1.0 (14) 0.6 (8) 365 (5.7) 14 (0.2) 1.2 

20 163 (2.6) 102 (1.6) 1.4 (20) 1.1 (16) 360 (5.6) 23 (0.4) 1.1 

25 158 (2.5) 83 (1.3) 1.5 (21) 1.3 (18) 355 (5.6) 28 (0.4) 1.1 

30 182 (2.9) 124 (1.9) 1.3 (18) 1.0 (14) 361 (5.7) 23 (0.4) 1.2 
" Counts degrec^^ (lO^'' Lq kpc"'^). 
Degrees (kpc). 

Counts degree~^(10~^ Lq kpc"-^). 
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Fig. 1. — A tangent-plane projection of the PAndAS fields around M33. The central field, 
m33c (black), is from the archive (the data is from 2003B and 2004B). All other fields with 
prefix m (red) were observed in 2008B. Fields with prefix nb (green) and tb (light blue) were 



observed in 2009B and 2010B, respectively. The dark blue solid e. 



(73 X 45 arcminutes) at which ^ 25 mag arcsec ^ (INilson 



lipse marks the diameter 



19731 ). The two perpendicular 



lines s 


low t 


Nilson 


1973) 



dashed-line circles mark radii of r = 1, 2, 3 and 3.75 degrees (14.1, 28.2, 42.4 and 53.0 kpc, 
respectively). We assurn e a di stance modulus of (m — Mq) = 24.54 ± 0.06 (809±24 kpc; 



McConnachie et al. 



2004 . 



20051 ). The three straight black dashed lines each represent one 



line of equivalent Galactic latitude (6 = -35.3, -31.3 and -27.3). 



- 46 - 




Fig. 2. — Color- magnitude (Hess) diagrams of the different annuli shown in Figured! Bins 
are 0.025 x 0.025 mag, and are shown with a logarithmic scaling in number counts of stars. 
Contamination due to the foreground MW halo and disk stars is estimated with the regions 
defined by 0.1 < - i)o < 0.6, 19 < io < 22, and 1.5 < - i)o < 3, 17 < zq < 20, 
respectively (shown as the boxes in each panel). The isochrones correspond to [Fe/H] = -1.0 
and -2.5 dex for a 12 Gyr, [a;/H]=0.0 stellar population at the distance of M33, and have 
magnitude limits of 21.0 < io < 24.0. The annulus between 3 < r < 3.75 degrees was used 
to determine the levels of foreground contamination. The bright clump at io ~ 25, < io < 
1 is mainly composed of misclassified background galaxies (with a very small number of M33 
horizontal-branch/red-clump stars) . 



^ (degrees) 



3 . — Density contours of candidate RGB stars similar to Figure 13 in 



McConnachie et al. 



20101) but updated using data from 2010B for frames tb62-tb66 (see Figured]). The grey 
contour is la above the background, corresponding to an estimated surface brightness hmit 
of fiv = 33.0 mag arcsec"^. We exclude regions within this contour for our estimate of the 
stellar halo. The black contours correspond to 2, 5, 8 and 12cr above the background {fiy ~ 
32.5, 31.7, 31.2, and 30.6 mag arcsec"^, respectively). The feature at (^, r] = (-3.5, 3.6) is 
Andromeda II. 
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Fig. 4. — The radial profiles of extended disk substructure (red triangles) and non- 
substructure (blue crosses) regions, identified by the grey \a contours in Figure El normalized 
using the total annulus area. The total radial profile is shown by the black squares. We 
only show the profiles beginning at ~0.85, within which the extended disk substructure com- 
pletely dominates. Each bin location is fixed for all three components, and all bins have a 
fixed width of 0.1 degrees (shown by the horizontal error bars). The vertical error bars show 
^Jnj area. 




Fig. 5. — Smoothed maps of the spatial distribution maps for the candidate galaxies, MW 
disk, MW halo and M33 RGB stars. See Section H] for details. Clearly visible in the RGB 
map (-2.5 < [Fe/H] < -1.0 dex) are the M33 extended disk substructure and Andromeda 
II in the NW. The ellipse, two perpendicular lines and circle are as in Figure [H To ensure 
the most effective colour range to show features (or lack thereof), zeropoints were set at log 
counts = 2 for the galaxy, disk and halo plots, and log counts = 4 for the RGB plot. 
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Fig. 6. — The left- and right-hand columns show the azimuthal and Galactic latitudinal 
distributions, respectively, of the density (counts degree"^) variations for each region in 
Figure [2l The data within the 3 < r < 3.75 degree annulus, having excised the area 
associated with the extended disk substructure, is shown. For the azimuthal distributions, 
0, 90, ±180, and -90 degrees correspond to east, north, west and south, respectively. M3rs 
centre is at approximately 135 degrees in this orientation (as indicated by the dashed line). 
The errors in all panels correspond to the values of -y/n/area. 
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Fig. 7. — The background-uncorrected (upper dashed curved hue and points) and the 
background-corrected (lower sohd curved hne) radial profiles of RGB candidate stars. The 
horizontal dashed line indicates the background level, S^g. Two radial density (vertical) 
error bars are shown: the smaller set is calculated using -y/n/area as the error in each bin. 
The bin size was allowed to vary until the required signal-to-noise ratio of 25 was reached. 
Horizontal "error" bars show the width of the bin. The vertical dashed regions indicate the 
radius within which we do not have any data because we excise the area dominated by the 
disk and extended disk substructure surrounding the disk. The larger vertical error bars 
show the variation due to residual substructure (see Section H73l for details). 
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Fig. 8. — The cumulative luminosity as a function of radius. Horizontal "error" bars show the 
width of the bin, and are the same as those shown in Figure [71 The error on the luminosity 
is calculated by combining the Poisson errors for the RGB candidate star counts with the 
uncertainty of the background. This figure highlights the increase in our estimate of the 
luminosity uncertainty as we increase the area we consider. 
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Fig. 9. — The azimuthal distribution of RGB candidate stars for the annuh shown in the 
top left of each plot. 0, 90, ±180, and -90 degrees correspond to east, north, west and south, 
respectively. M31's centre is at approximately 135 degrees in this orientation (as indicated 
by the dashed line). The top, middle and bottom panels show the data with r < 3, 1 < r < 
2, and 2 < r < 3 degrees, with the regions associated with the extended disk substructure 
excised in all panels. The left- and right-hand arrows correspond approximately to the SE 
and NW tips of the S-shaped warp of the extended disk substructure. 



0.0 0.5 1.0 1.5 2.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 2.0 
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Fig. 10. — Similar to Figure Ej but here showing the CMD for the region with r < 3 degrees. 
The left-hand and middle panels show the CMD after excising the extended disk substructure 
areas in Figure |3] (which effectively imposes a minimum radius of r = 0.88 degrees). The 
RGB that we aim to detect is so faint that it is barely visible on the left hand plot. We 
overlay a [Fe/H] = -2 dex isochrone on the middle plot. The right-hand plot shows the CMD 
of the extended disk substructure areas for comparison. 
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Fig. 11. — Background-uncorrected profiles for the quadrants split by major and minor axes, 
e.g., as shown in Figure [31 




Fig. 12. — The CMDs for the quadrants used in Figure fTTl 



-57- 




Fig. 13. — The halo luminosity as a fraction of the host galaxy luminosity against the total 
host galaxy mass (dark plus luminous) for the MW, M31, M33 and NGC 2403. Data for the 
MW, M31 and NGC 2403 are from the literature (see Secti ons [211 and 151 respec tively, 



for details and references). M33's host galaxy mass is from iCorbelli fc Saluccil (120001 ). The 



ranges of M33's halo luminosity are our estimates from this paper. The observed range comes 
from strictly limiting the integration between the range of our actual data, whereas the 
extrapolated range comes from extrapolating inwardly to the centre of M33 and outwardly 



to M33's viria^ 



Purcell et al. 



radius. The lines represent the models for the intrahalo light fraction in 
(120071 ). We use Uefj = 1 in conjunction with the three values for shown 



in the plot. Ueff represents the effective number of satellites with mass Mgat = Mhostf^O', 
fd represents the total stellar mass fraction a satellite contributes to its host galaxy halo. 
No distinction is made between halo substructure and a smooth halo component for the 
host lu minosity estim ates in th e literature for the MW, M31, and NGC 2403 - similarly 
for the 



Purcell et al. 



models. 



McConnachie et al. 



(|2010[ ) estimate the luminosity of the 



extended disk substructure (EDS), Leds ~ O.OILmss so Leds ~ lO^L0.M33's extended 
disk substructure is in the halo region, and if it was to be included in the estimates shown 
in this figure it would raise each Lhaio/ Lhostgaiaxy estimate by 0.01. 



